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Abstract 



The final stage of Bose-Einstein condensation in a large volume occurs 
^ ' through coarsening — growth of individual correlated patches. We present an- 

. alytical arguments and numerical evidence that in the momentum space this 

growth corresponds to strong turbulence of the particle number, with a turbu- 
lent cascade towards the infrared and the power law n{k) oc in d= 2,3 
■ spatial dimensions. As a corollary, we find that the correlation length grows 
linearly in time, with the speed of order of the speed of Bogoliubov's quasipar- 



ticles (phonons). We use these results to estimate the speed of EEC in atomic 



Q_i' vapors and in galactic dark matter. 



bep-ph/0201163 



1 Introduction 

The discovery of Bose-Einstein condensation (BEC) in alkali vapors had stim- 

ulated interest in a theoretical description of non-ideal Bose g research topic 

that dates back to Bogoliubov's work 0] on the equilibrium properties. In current 
experiments, the gas is confined in a small magnetic trap and is strongly affected 
by the confining potential. However, theorists have also investigated the kinetics of 
BEC in a uniform (on average) gas, contained in a macroscopic volume 0-11^ • In 



addition to the intrinsic interest that this problem has for basic science, it is also of 
potential importance for astrophysics and cosmology in connection with the problem 
of galactic dark matter [|T5|-[|19|. 

The earliest theoretical work we are aware of on this subject is the paper by 
Kraichnan ^ . Although we will confirm below several qualitative results of that work 
(such as concentration of particles in a progressively smaller momentum range and a 
turbulent cascade at larger momenta), the quantitative details differ. We attribute 
the difference to the imprecise treatment of phonons in ref. (noted by Kraichnan 
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himself). Although the amount of theoretical work on BEC in a macroscopic volume 
has grown rapidly in recent years an adequate understanding of the phonon 

dynamics has remained elusive. In the present paper, we attempt to improve the 
situation, through a combination of analytical estimates and numerical results. 

We consider an isolated nonrelativistic Bose gas, sufficiently dilute for the inter- 
action between the particles (atoms in the laboratory, or dark matter particles in the 
outer space) to be characterized by a single parameter — the scattering length. Such 
a gas can be described by an effective Hamiltonian: 

where and ip'^ are canonically conjugate (field) variables, m is the particle mass, and 
Qi is the coupling that encodes the value of the scattering length.|^ In a macroscopic 
volume, a negative (74 (i.e. an attractive interaction) leads to clumping, the process 
that had been investigated numerically in ref. [|1^]. In the present paper we consider 
the case of a positive (a repulsive interaction). In the concluding section, we will 
also discuss briefiy the case when the coupling is switched from positive to negative. 



as in a recent experiment ||2^. We assume that the volume is macroscopic, and the 
gas is on average uniform. 

A fundamental role in the kinetics of a Bose gas is played by the length scale 
corresponding to the wavenumber /cbo, 

^^^Bo = ^9imn^^e , (2) 

which we will call the Bogoliubov wavenumber. Here riave = (V'^V') is the average 
density of the gas. Particles with momenta p ^ ^^bo can be considered as weakly 
interacting, while those with smaller momenta in general cannot. For p <C hk-^o, 
elementary excitations of a fully condensed gas have a linear dispersion law; these are 
Bogoliubov's phonons 

In addition to the Bogoliubov length Abo = Stt/Zeboj which is time-independent, 
let us consider the time-dependent correlation length of the gas Ao(t). The corre- 
lation length is defined as the characteristic scale at which the two-point function 
('?/'''^(x)'?/'(0)) changes its behavior (e.g. from a power law to an exponential). The cor- 
responding wavenumber is koit) = 27r/Ao(t). BEC corresponds to the growth of the 
correlation length with time, and we now summarize some previous results concerning 
this growth §-[0. 



If the initial value of the correlation length is small (we say that the initial sate is 
disordered) , 

Ao(0) < Abo , (3) 

BEC proceeds in two stages. At the first stage, Ao(t) ^ Abo- If one assumes that 
the main role is played by interactions of particles with wavenumbers larger than or 

* In the context of dark matter, using (|^) implies neglect of gravitational interactions, which in 
general is only possible up to some spatial scale. We discuss this question further in the conclusion. 



2 



of order ko{t), one can use a quantum Boltzmann equation, because at this stage all 



such particles are weakly interacting. The solution to the Boltzmann equation |]TT 
shows that the main feature of the spectrum is a turbulent front that propagates to 
smaller wavenumbers. The position of the front can be identified with kQ{t). 

When the correlation length exceeds Abo, the particles begin to interact strongly 
with their local environment, and the standard Boltzmann equation cannot be used. 
At this second stage the growth of the correlation length continues through coarsening, 



i.e. a process of alignment of the field on progressively larger spatial scales 
characterized by a power-law dependence of Aq on time (at large times): 

Ao(t) = const X . (4) 

We will return to discussion of the scaling exponent a shortly. Also, as we will see 
below, the particle number contained in modes with k < /cbo is in fact concentrated in 
the ever-shrinking region k ~ fco(^)- So, at this stage the correlation length Aq gives 
the typical size of correlated patches (i.e. regions over which the phase of the field is 
approximately constant). 

If the initial value of the correlation length is already large (the initial state is 
partially ordered), 

Ao(0) > Abo , (5) 

there is no room for a "preparatory" stage that can be described by a Boltzmann 
equation for the particles: one has to use a different description from the beginning. 
Here we consider only the late-time coarsening stage, which is present in either of the 
two cases. 

How does one go about determining the scaling exponent a? The fastest the 
correlation length can possibly grow is at the speed of sound (cf. ref. 0), which in 
our case is the speed of Bogoliubov's phonons 

Vs = {gAn^.^e/mY^'^ ■ (6) 

If one replaces (|^) with a theory of non-interacting phonons, as it is done in the 
"phase-only" model of ref. |T^, the only scaling law one can possibly discover is 
Ao(t) ~ Vst, corresponding to a = 1. Although this value of a is consistent with 
the results of numerical simulations ||12| of the full theory (|l]), the questions remain 
if the approximation of non-interacting phonons is a good one, and why their free 
propagation should have anything to do with the growth of the correlation length. 

In the present paper, we attempt to answer these questions by taking a look at 
what happens in the momentum space, i.e. at the evolution of the power spectrum 
of the field. We find that a consistent description of the power spectrum is obtained 
by assuming that the longest-wave phonons, far from being non-interacting, are in 
the opposite, maximally nonlinear, regime. This assumption leads to a prediction 
of an inverse (i.e. directed towards the infrared) turbulent cascade of the particle 
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number, with a definite value of tlie scaling exponent. This prediction is confirmed 
by numerical simulations in two and three spatial dimensions. 

We stress that this strong turbulence, characterized by strong quasiparticle inter- 
actions, is completely different from weak turbulence found in a variety of physical 
systems on the basis of a Boltzmann equations for weakly interacting quasiparticles 



21| . Likewise, it is different from the turbulent solution obtained numerically in ref. 

IT I for the Boltzmann stage of the evolution of a Rose gas. 

Because phonons at wavenumbers /c ~ I/Aq are maximally nonlinear, the rate at 
which the particle number is transferred to still larger scales is of order of the phonon 
frequency uj = Vgk. The time at which the order is established at a scale k is therefore 
t ~ l/w, which leads to 

Ao(t) ~ Vst . (7) 

Thus, the linear scaling law for the correlation length obtains, and the speed of EEC 
is of order f^, even though phonons are by no means free-propagating. 

The paper is organized as follows. In Sect. 2 we present analytical considerations, 
and in Sect. 3 the numerical results. Concluding Sect. 4 contains estimates of 
condensation times for trapped atomic gases and galactic dark matter. 



2 Acoustic turbulence 

For a fully ordered Bose gas (i.e. a Bose-Einstein condensate plus small fluctuations), 
the spectrum of quasiparticles had been found by Bogoliubov for an arbitrary 
pairwise potential. For the local interaction contained in (|I]) the dispersion law, oj{k), 
is given by 

uj (k) = k + -— . (8) 

m 4m^ 

Hence, the wavenumber k-Qo, eq. (|^), separates those k for which the spectrum is 
significantly modified by the interaction with the condensate from those for which 
quasiparticles approximately coincide with the ordinary particles of the gas. 

For a partially ordered gas, with the typical size of correlated domains (the 
correlation length) equal to some Aq, eq. (|]) can be applied to excitations with 
k ^ ko = 271 /Xq, because these can be viewed as propagating on an ordered back- 
ground (we will make the latter statement more precise shortly). In this section, we 
will present analytical arguments in favor of a certain form of the power spectrum of 
ip in the low-momentum range 

/CO < A; < ^Bo • (9) 

Recall that we consider the stage when ko <^ kso- For wavenumbers in the range 
(§), the quasiparticle dispersion law u){k) is nearly linear in k, so we refer to the 
excitations as phonons. 
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In perturbation theory, one writes ip = y/n exp{i6 / h) and defines the zeroth-order 
approximation by replacing n with nave in the first term in (|l|). This gives 

H^Ho = ^{Vef + ^g,n'. (10) 

From the kinetic term in the Lagrangian, ihijj^dtip, we find that —n is the momentum 
conjugate to 6. We now see that (|T0|) is a theory of noninteracting waves with the 
linear dispersion law 

ij\k) = sfh^k' . (11) 



This theory coincides with the "phase-only" model considered in ref. |T^, except 
that we do not introduce any coupling to a heat bath. Eq. ( pA]) is, of course, the 
low-momentum limit of (|^). 

Next, allow n to deviate from nave in the gradient term and consider small 6n = 
n — riave- To the second order in 6n, the Hamiltonian becomes H ^ Hq + Hi, where 

„,^^^V0)'^p-^. (12) 

2m 8m nave 

The second term results in a correction to the dispersion law, restoring the full eq. 
(§). (A similar way of deriving Bogoliubov's spectrum was used in ref. 0.) The first 
term is an interaction among quasiparticles. Using the decomposition of 6 and Sn in 
terms of quasiparticle creation and annihilation operators. 



^ - ^i:/||(«.+«ye-. (13) 

where V is the total volume, we find that the matrix element of the interaction, 
U(ki, k2, ks), has the conventional phonon form at low momenta: 

UikiMMocihhhf' . (15) 

We can now understand better the sense in which we can talk about quasiparticles 
in a partially ordered state. First, note that in addition to perturbative excitations 
there are nonperturbative ones — vortices in two dimensions and vortex lines in three. 
These cannot be described in terms of quasiparticles because n vanishes at the core 
of a defect, making it impossible to even define 6 there. We expect that the typical 
distance between the defects is of order of the correlation length Aq (cf. ref. |TO[)- 
Because perturbative excitations scatter strongly on the defects, the quasiparticle 
concept breaks down for k < ko: such excitations do not have any room to propagate. 
Of course, excitations with k ^ k^ also scatter on the defects, but only on scales of 
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order Aq. On shorter scales, they can be treated as if the defects were absent. In 
particular, if the interaction among quasiparticles themselves leads to a mean-free 
path shorter than Aq, it is that interaction (rather than scattering off the defects) 
that determines the quasiparticle lifetime. We will now argue that such is indeed the 
case. 

Indeed, it is precisely for acoustic (close to linear) dispersion laws that the con- 
ditions under which quasiparticle interactions can be treated perturbatively become 
quite subtle. The reason has been explained in ref. and is the following. For 
an exactly linear dispersion law, wave packets do not disperse, so if two excitations 
are produced with nearly parallel momenta, they will stay together for a long time, 
breaking the usual condition of applicability of perturbation theory. In our case, the 
dispersion law uj^k), given by eq. (||), is not exactly linear but curves upward from 
a straight line. As a result, the process in which two quasiparticles merge into one 
is kinematically allowed but, because the deviation of the dispersion from linear is 
small at small k, the initial momenta have to be almost parallel. Thus, quasiparticles 
are forced into the kinematic arrangement in which the interaction between them 
can be large. In particular, a weak turbulence spectrum, obtained as a solution to a 
quantum Boltzmann equation for phonons (with an interaction of the form (|l^)), is 
applicable only at sufficiently large k |^3[. Therefore, in the small k region (^), we 
will be searching for strong turbulence, for which the lifetime of a quasiparticle is of 
order of its inverse frequency, and the Boltzmann equation cannot be applied. 

To see further the importance of the scattering process in which two quasiparticles 
merge into one, let us visualize BEG as going from the situation shown in Fig. |I]a to 
that shown in Fig. ||b. The field ip is represented by arrows on a grid, with the length 
of an arrow corresponding to and the angle of rotation to aigip oc 6. The grid is 
made one- dimensional for illustrative purposes only; in fact, we do not consider d = 1 
in this paper. We have chosen these particular configurations in such a way that the 
lengths of all arrows are about the same, and the magnitudes of angular differences 
between neighboring arrows are also the same, for all sites in Figs. |l|a and |l|b. Thus, 
both the gradient and potential energies are about the same for the two figures, i.e. 
the energy is conserved. Nevertheless, the state of Fig. |^a has no long-range order, 
while the state of Fig. clearly has some. We can also say that the first state has 
many low-momentum quasiparticles, while the second state has relatively few high- 
momentum ones. The transition from the first state to the second should therefore 
occur through merging of quasiparticles, leading to a decrease in their number and 
an increase in the average quasiparticle energy. 

We now see that the process of BEC is largely analogous to turbulence in an 
incompressible fiuid. Instead of turbulent eddies we now have quasiparticles, and 
instead of a larger eddy splitting into smaller ones we now have two quasiparticles 
merging into one (of course the inverse process exists as well but, as we have seen, 
for BEC to occur the merging of quasiparticles should statistically prevail over the 
splitting). Moreover, we have seen that the interaction among quasiparticles in this 
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t/t\t/t\t/t\ 

(h) 



Figure 1: Difference between uncondensed (a) and condensed (b) states with equal 
energies and equal particle numbers. 

process has every chance to be strong, so that the corresponding mean-free path can 
be short — shorter than the correlation length Xq. If that is indeed so, this process is 
the main interaction effect for the entire momentum range @, so this range becomes 
analogous to the inertial range in an incompressible fluid (cf. ref. We can 

therefore try to deduce the power spectrum of ip by an argument a la Kolmogorov 
and Obukhov (as reviewed for instance in the Landau-Lifshitz textbook [Q). 



Unlike the case of an incompressible fluid, where only the energy is conserved, in 
our case there are two conserved quantities: the energy and the number of particles (by 
particles we mean here the ordinary particles of the gas, as opposed to quasiparticles, 
whose number is not conserved). The process of BEC can be viewed either as particle 
number accumulating in low-momentum modes or, as we have seen above, as energy 
accumulating in high-momentum modes. So, we have to decide with respect to which 
quantity the turbulent cascade will be established. A posteriori checks show that the 
only consistent choice is an inverse (i.e. directed towards the infrared) cascade of the 
particle number. 

We therefore proceed to deflne the particle number contained in all scales shorter 
than a certain scale A = 27r/k, with k in the inertial range (^: 

N{k,t)= J2 (4'^k'(t)), (16) 

|k'|>fc 

where the angular brackets denote averaging over the quantum state of the system. 
A turbulent cascade corresponds to a constant (fc-independent) flux of the particle 
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number in momentum space. The flux can be estimated by multiplying A^(A;) by 
the rate at which quasiparticle interactions transport the particle number to larger 
scales. This rate is simply the inverse quasiparticle lifetime; for strong turbulence it 
is of order of the quasiparticle frequency uj{k) ^ Vgk. Thus, we obtain 

kN{k,t) = A{t) , (17) 



where A is a /c-independent function of time. Eq. (|T7|) corresponds to the following 
power spectrum of ip 



2^ m 



where A is proportional to A, and ci = 2 or 3 is the number of spatial dimensions. 

We can now consider some consistency checks. According to the previous dis- 
cussion, we expect that the particle number, given by the integral over k of P{k) 
saturates in the infrared, while the gradient energy, given by the integral of k'^P{k) 
saturates in the ultraviolet. Both conditions indeed hold for the power spectrum ([T8|) . 

The integrals in question are cut-off at the boundaries of the range (^). In partic- 
ular, the integral corresponding to the total particle number is cut-off at /c ~ /cq: 

Ntot - , (19) 

where V is the total volume. The time dependence of k^ can be found by pushing our 
estimate for the quasiparticle lifetime T{k) ~ 1/vsk to the limit of its applicability, 
i.e. to A; ~ /cq- Then, r(/c) can be interpreted as the timescale at which the turbulent 
spectrum reaches the wavenumber k, or equivalently the size of correlated domains 
extends to spatial scale A = 27i/k. The estimate t(A;o) ~ l/'^^s^o corresponds to 

Ao(t) ~ v,t , (20) 

i.e. a linear growth of the correlation length with time. We stress that this growth 
is due to a strong interaction among quasiparticles and cannot be described by the 
phase-only model (0), which contains no such interaction at all. 

For comparison, in the ring-model approximation for d = 3 Kraichnan finds 
P{k) oc k~^^^^ (which is curiously close to our scaling ([l8|) ), and Ao(t) oc (Note 



that notations differ: our k^ corresponds to Kraichnan's k^,, and our hk'^^/2m is twice 
his Uq.) 

Because the total number of particles is time-independent, the scaling ko oc 1/t 



implies, through (0), that 

, , , const , , 

m = . (21) 



Thus, the spectrum propagates to smaller k according to (poD, while the magnitude 
of the spectrum at a given k decreases according to (pT]). 
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We wish to stress again that the power spectrum (|T^) describes the distribution 
of ordinary particles of the gas, as opposed to quasiparticles. This distinction should 
be kept in mind when making comparisons with the literature. For instance, weak 
turbulence is often described by a distribution of quasiparticles ||2l| , the same ones for 
which the Boltzmann equation is written. In our case, a transition to quasiparticle 
spectra may be difficult: due to the strong interactions, the statistics of quasiparticles 
is likely to be non-Gaussian and therefore not amenable to a simple description in 
terms of occupation numbers. 



3 Numerical results 

The results of the previous section are all based on one central assumption, namely, 
that a turbulent cascade, described by a constant flux of particles in the /c-space, is 
indeed established. As well known in other contexts, a turbulent cascade can only 
occur if the interaction (in our case, among quasiparticles) is in a certain sense lo- 
cal in the /c-space. In our case, that means that the main role should be played by 
interactions among quasiparticles whose momenta are of the same order of magni- 
tude, rather than, say, by processes in which a quasiparticle will jump directly to a 
momentum of order /cbo- 

Another ingredient of our theory is the use of the maximally nonlinear estimate 
for the quasiparticle lifetime. Although we tried to motivate this estimate in the 
previous section, it should still be properly considered as an assumption. 

For cases of weak turbulence, the assumption of locality can be verified, and the 
quasiparticle lifetime can be found by using the solution to the corresponding Boltz- 



mann equation We are not aware of any such analytical method for strong 

turbulence. So, we have resorted to numerical simulations. The test of the assump- 
tions will be indirect — through the form of the power spectrum. 

Operationally, the numerical approach is quite straightforward. One assumes that 
the pile-up of particles at low momenta leads to large occupation numbers in low- 
momentum modes, so that the field ip can be regarded as classical and described 
by a classical equation of motion (which for the Hamiltonian (|T]) is known as the 
Gross-Pitaevsky equation). The classical equation can be applied to the coarsening 
stage of the evolution |[T^ (as well as to the Boltzmann one |g]), as long as the 
assumption of classicality is good. This approach was used in the previous numerical 
study of BEG [0, but the form of the power spectrum was not addressed there. (For 



attractive interactions, the classical approximation was used in refs. |I3| , p!^.) Here 
we concentrate on the power spectrum of ip. 

The condition for the occupation numbers to be large, though, is not entirely 
trivial. For wavenumbers in the range (^, it can be derived if we use the already 
obtained power spectrum ( |TB| ) (anticipating that this form of the spectrum will be 
confirmed numerically). Defining the particle occupation numbers as = (V'k'^k) 
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and using (|19]), we find 

^avcfco(t) 

• (22) 

Thus, at any given k the occupation number actually decreases with time. Still, for 
all k from the range 



^a vefeo(t) _ fco(t) 

where e = k^^/n^ve- So, if the interaction parameter e is small, the condition rik ^ 1 
can hold even for an already small ratio /cq/^Bo- 

Furthermore, the use of the classical approximation implies that ^ 1 for all 
"important" momenta. That is certainly not so for a thermal state, in which both 
the particle number and the energy are concentrated in modes with nk ~ 1. For this 
reason, we will consider nonthermal initial conditions. In addition, we will have to 
make sure that the high-momentum modes do not get overly populated in the course 
of the evolution. This latter requirement is rather nontrivial, since the gradient energy 
contained in the range (H) is, for the power spectrum ([T8|), of order 



Z*"^" dkeP{k) ~ A(t)fcBo , (24) 

Jfco 



and thus decreases with time. So, unlike the particle number, the gradient energy is 
not contained in the low momentum modes but is being transferred to modes with 
k > k-Qo- In this way, the low-momentum part of the gas gets "cooled down", while 
the high-momentum part is "heated up". In a quantum system, the flow of energy 
to high momenta leads to thermalization, but in a classical system on a lattice, the 
energy eventually becomes dominated by the cutoff modes. When that happens, the 
simulation breaks down. Fortunately, when we plot the power spectra, we will be 
able to monitor their high-momentum content as well. 

After these remarks on the applicability of the classical approximation, we proceed 
to solve the classical equation of motion for a dilute Bose gas (the Gross-Pitaevsky 
equation), which after some rescalings takes the form 

= -V'ij + . (25) 

at 

The rescalings that lead from (|l|) to (^) are as follows. A factor of is absorbed 
into ijj, a factor of h into the time derivative, and a factor of h/\^2m into the spatial 
derivatives. We hope that using the same letters for the rescaled quantities as for 
the original ones will cause no confusion, as only the rescaled quantities are used 
throughout this section. The rescalings correspond formally to setting h = 1, = 1, 
and m = 1/2. In particular, the expression for the Bogoliubov wavenumber becomes 

^Bo = A/2nave ■ (26) 
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The power spectrum of a classical field is defined by averaging over direction of the 
wavevector: 

Pik) = 7^ / dnM' , (27) 

"tot J 

where f2tot = 27r in c? = 2, and Qtot = 47r in (i = 3. It satisfies the equality 

-d / rd 



n,,e= i2n)-'' J rkP{k) . (28) 

Eq. ( p5| ) was solved in two- and three-dimensional boxes will equal sides of length 
L and periodic boundary conditions. We used an operator-splitting algorithm, based 
on alternation of two different types of updates corresponding to the individual terms 
on the right-hand side of (]25|). The updates were alternated so as to maintain the 
second-order accuracy in time. The action of was computed in the momentum 
space, using the fast Fourier transform. The algorithm conserves the number of 
particles exactly. Energy nonconservation, for the times shown in the figures, was 
about 0.2% in = 3 and about 2% in = 2, relative to the value of energy at t = 0. 

The initial condition was taken in the form 



V'k(t = 0) = ^P(A;,0)exp(^7k), (29) 
with random uniformly distributed phases 7k and a nonthermal power spectrum 

P{k, 0) = B exp(-A;V/i2) . (30) 

For both runs shown in the figures, we used = tt, while the values of B were such 
that the average densities equaled riave = 31.4 for d = 2, and riave = 27.8 for d = 3, 
corresponding to k^o = 7.93 and fcso = 7.46. The parameter /i can be viewed as 
the wavenumber corresponding to the initial correlation length Ao(0): /i = vr means 
Ao(0) = 2. Because fi < /cbo ioi both runs, the gas was in the coarsening regime from 
the beginning. 

In d = 2, a Bose gas cannot support long-range order at any finite temperature 
However, a classical simulation of an isolated system described by eq. (P^D, with 



a finite total energy and a large (ideally, infinite) ultraviolet cutoff, corresponds to 
an effectively zero temperature. In particular, at a finite temperature there would 
be a finite equilibrium density of vortices (or vortex-antivortex pairs). In a classical 
simulation, the number of vortices and antivortices decreases with time. 

We show the evolution of the power spectrum for c? = 3 in Fig. ^ and for d = 2 in 
Fig. ^. We see that in either case a region of power law forms, well approximated by 
eq. (0). According to our previous discussion, the power law should extend in the 
ultraviolet to wavenumbers of order /cbo- Numerically, we find a "knee" in the power 
spectrum at ~ 2fcBo- 



It is interesting to compare our results to those obtained in ref. for an exter- 
nally driven d = 2 Gross-Pitaevsky equation. In that case, the total number of parti- 
cles is not conserved, so the parameters A;bo and Vs are time-dependent. Nevertheless, 
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Figure 2: Power spectrum of ^ in = 3 at three different moments of time, obtained 
on a 64^ lattice witli side lengtli L — 9. The straiglit line is a power law. 
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Figure 3: Power spectrum of in d = 2 at three different moments of time, obtained 
on a 256^ lattice with side length L — 9. The straight line is a k"^ power law. 



12 



if the change of these parameters is sufficiently slow, we can still expect formation of 
an inverse cascade via quasiparticle scattering.^ The authors of ||2^ obtain a power- 
law spectrum oc A;~^, i.e. an index different from ours, but they also observe that 
the spectrum (including the index) is anisotropic. We believe the anisotropy can be 
explained by the presence of some remaining vortex-antivortex pairs, perhaps just 
one such pair in the entire integration volume. Such well-separated vortices are easily 
pinned on a lattice, so a single pair can survive for a long time, causing an effective 
anisotropy. Then, it makes sense to average the power spectrum over the direction of 
the wavenumber, as we do in our definition (PTj). In the presence of anisotropy, the 
main contribution to the average will come from directions in which fiuctuations are 



the largest. It is precisely in such directions that the spectrum obtained in ref. p6 
is very similar to ours, with a relatively fiat region at the smallest k, and a power law 
nk oc at larger k. 

The spatial sizes of our lattices were not sufficient to perform quantitative checks 
of the temporal scaling laws (EO) and (pT]). While for the first of these one may be 



content with the results of ref. it would certainly be interesting to directly test 



(pip on larger lattices. Qualitatively, the predicted by (|2TD decrease in the magnitude 
of the power spectrum in the inertial range is clearly seen in both Fig. |^ and Fig. ^ 



4 Discussion 

We now apply our results to two cases: (i) BEC in atomic vapors and (ii) formation 
of a Bose-Einstein condensate from galactic dark matter. In both cases, we estimate 
the speed of BEC as being of order of the sound velocity @. We can then estimate 
the time scale required for BEC to complete on physically interesting spatial scales. 

Using the estimate 4 x 10^^ cm~^, obtained in ref. for the number density of 
sodium atoms, as nave in eq. and noting that our parameter is identical to U 
used in that paper, we find 5'4nave = 520 nK, and 

Vs = 1.4 cm/s . (31) 

Thus, for a gas cloud a few micron across, we estimate that the coarsening time is 
shorter than 1 ms. The duration of evaporative cooling in ref. was 7 s, so the 
coarsening time was not a limiting factor in that experiment. In the experiment with 
rubidium [|^, density was about two orders of magnitude smaller, which reduces Vg 
by about an order of magnitude compared to ( [HT| ) {g^ and m are also different, of 
course, but by smaller factors). Still, the coarsening time remains much smaller than 
the time scale of evaporative cooling. 

It is not immediately clear what the rate of change of A'tot should be compared to. However, 
for the case called moderate pumping in ref. this rate is much smaller than even the smallest 
energy scale present in the Hamiltonian (|l|), h fcj„;,j/2rn, where fcmin is the minimal nonzero lattice 
wavenumber. It seems plausible that in this case the quasiparticle spectrum and the scattering rates 
are not affected too much by the time-dependence of iVtot- 
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On the other hand, the relaxation time corresponding to the first, Boltzmann, 
stage of BEG in these gases is quite large. For a gas of hthium-7 atoms (which have 
an attractive interaction at low densities), with the parameters of the experiment of 
ref. 0, this time was estimated in ref. |]13[ to be of order 1 s. Similar estimates 
can be obtained for gases with repulsive interactions. First, form the dimensionless 
parameter of nonlinearity 

ml 

where a is the scattering length, and T is the temperature in energy units. For 
rubidium-87, the authors of ref. [0 quote a ~ 10^^ cm, n = 2.5 x 10^^ cm~^, and 
T = 170 nK. With these values of the parameters, C,t ~ 10~^. In the Boltzmann 
approximation, the relaxation rate is proportional to a^. Hence, we can roughly 
estimate the relaxation time as 



h 



tr ~ 7^ ~ 0.4 s . (33) 



The authors of |T| noted that after reaching the condensation temperature the con- 
densate fraction grew for 2 seconds. Given how crude the estimate ( p3| ) is, it seems 
plausible that this delay can be attributed to the Boltzmann relaxation. 

The distinction between disordered and partially ordered initial states is important 
also for a Bose gas with an attractive two-particle interaction (a negative a). The 
wavenumber /cbo is now defined (in units with ^ = 1) by 

^Bo = 4|fi'4|mraavc = 167r|a|nave (34) 

(the last equality is specific for three spatial dimensions). The quadratic dependence 
of the rate of collapse (fragmentation into individual drops) on the scattering length 
13| applies only when the initial state is disordered, i.e. the correlation length Aq 



is smaller than Abo- In the opposite limit of a large correlation length, collapse in a 
sufficiently dense system will occur by linear instability. Indeed, the dispersion law 
for excitations with ^ fco is now 

uj\k) = ^,{k'- kLk') , (35) 

which corresponds to an instability for all k < /cbo- The instability is the strongest 
for k = k-Qo/ v^, where its rate is 

l^^lmax = /^Bo/4"^ ■ (36) 

A large initial correlation length can be achieved, for example, if the scattering length 
was originally positive, but after some time was switched to a negative value, as in a 
recent experiment . If we imagine that the gas is confined in a sphere of radius R 
with Dirichlet boundary conditions for ip, and we gradually increase /cbo, by increasing 
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\a\ or riave, the instability will first appear when k-Qo = tt/R. Then, after kso passes a 
transitional region with k^o ~ 7r/i?, the time scale of collapse will, as seen from (|36|) , 
become inversely proportional to the scattering length and to the density of atoms. 
In our estimation, these scaling laws — especially the former — are well born out in the 
experimental results of ref. PO . 



For Bose gases considered in cosmological and astrophysical applications, one 
should take into account gravitational effects. If the particles are quanta of a real 
scalar field with a potential 

1^(0) = ^^^0' + , (37) 
then, in the absence of gravity, (0) is the correct nonrelativistic limit, provided 

This limit is obtained by writing 

= Ue-'""' + T^^e™*) (39) 



'2m 

and assuming that i/j changes with a frequency scale much smaller than m. Eq. (^) 
for the sound velocity becomes 



1/2 

(40) 



where p = mnave is the mass density. 

Gravitational interaction (in Newtonian gravity) modifies the quasiparticle dis- 
persion law, so that now 

u^k) = {k' + kl^e - kj) , (41) 

where kj = 167rGm^nave is the Jeans scale of the gas. (For a discussion of various 
length scales in a self-gravitating Bose gas, and various types of gravitational insta- 
bilities, see ref. H^-) In a partially ordered Bose gas, eq. (0), similarly to eq. (§), 
applies only at scales smaller than the correlation length Xo{t). It follows that once 
Ao exceeds a critical value Amax, 

kj 4Gm^ ' ^ ' 

the gas becomes gravitationally unstable via a linear (Jeans) instability and fragments 
into smaller clumps. Similar expressions were obtained in [^, [|T^-[|T^ for the radius 
of a stable gravitationally bound object. 
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According to the general theory presented above, BEC at scale Amax will occur 
through coarsening only if Amax is larger than the Bogoliubov wavelength Abo, or 
equivalently /cbo > kj. This corresponds to 

-25 ( m Y 




10-^— J . (43) 



In the last estimate, we have used p = O.O2M0/pc'^ for the mass density at a 
galaxy core. The time scale r required for coarsening to extend to scale Amax can be 
estimated as 

Thus, regardless of the coupling, r is much smaller than the current age of the universe 
(although much larger that the age at the time of galaxy formation). 

On the other hand, the condition for the first, Boltzmann, stage of BEC to com- 
plete, obtained in |]16|, |l^] is, for the same value of p, 



If this condition is satisfied, completion of BEC on scale Amax is guaranteed by (|4^). 
So, for the gas to be at least linearly stable on the scale of the galaxy core, Amax 
should exceed the core diameter 2rc. In view of (|4^) , this requires 



Except for very small masses, ( ^6] ) is a stronger condition than (^5|) for all scales 
of relevance to galactic dark matter, Tc > 1 kpc. If satisfies (|^) but not (|46|) , 
then the coarsening eventually sets off a Jeans instability, and the core fragments into 
smaller objects (Bose stars). 

In conclusion, we have presented analytical arguments in favor of strong turbu- 
lence of the particle number occurring at the final, coarsening stage of Bose-Einstein 
condensation. We have found the form of the spectrum in the region of wavenum- 
bers corresponding to an inverse cascade of the particle number. We have confirmed 
both the existence and the scaling exponent of the power-law spectrum by numerical 
simulations of classical fields in two and three dimensions. According to our theory, 
quasiparticles (Bogoliubov's phonons) are in the maximally nonlinear regime, imply- 
ing that the correlation length grows linearly in time, with the speed of order of the 
speed of phonons. Finally, we have applied these findings to estimate the time scale 
of BEC in atomic vapors and in galactic dark matter. 
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